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Abstract. We present a compilation and study of DENIS, 2MASS, ISOGAL, MSX and IRAS 1-25/im photometry 
for a sam ple of 441 late- type stars in the inner Galaxy, which we previously searched for 86 GHz SiO maser 
emission jMessineo et alj 2002). The comparison of the DENIS and 2MASS J and Ks magnitudes shows that 
most of the SiO targets are indeed variable stars. The MSX colours and the IRAS [12] — [25] colour of our SiO 
targets are consistent with those of Mira type stars with dust silicate feature at 9.7 /xm feature in emission, 
indicating only a moderate mass-loss rate. 

Stars: variables: general - circumstellar matter - masers 



Il996|) . AGB stars thus permit a study of Galactic kine- 
matics, structure, and mass-distribution. 

To understand the Galactic structure and kinemat- 
ics it is important to combine the kinematic informa- 
tion and the stellar properties, e.g. luminosities, which 
can provide a distance estimate. Good photometry on in- 
frared point sources toward the inner Gala xy is now avail- 
able from larg e surveys su ch as DENIS tepcht ein et al, 
1994), 2 MASS llGutri et alf2 003b ISOGAL ( Omont et al, 
2003;: iSchuUer et al.ll200.C and MSX fegan et alJ Il999l: 
Price et al I2OOII) . Since the high extinction toward the in- 
ner Galaxy precludes studies at optical wavelengths, these 
infrared data permit a unique view of its stellar popula- 
tion. The combination of near- and mid-infrared photom- 
etry enable us to examine the nature of the stars, i.e. to 
derive their luminosities, mass-loss rates, and to discrimi- 
nate againts foreground stars. 

To improve the line-of-sight velocity statistics, we con- 
ducted 86 GHz i; = l,J = 2^1 SiO maser line 
observations of 441 late-type stars in the inner Galaxy 
(30° < I < -4° 161 < 1) with the IRAM 30-m telescope 
llMessineo et al.lE002l hereafter Paper I). This paper is 
part of a series devoted to characterise the properties, i.e. 
mass-loss rates and luminosities, of the 441 sources previ- 
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1. Introduction 

Stars of intermediate mass, 1 < Af* < 6 M©, enter a phase 
of intense mass loss when they reach the Asymptotic Giant 
Branch (AGB). As a consequence, they are surrounded 
■ with a dense envelope of dust and molecular gas. Due to 
I the low effective temperature and the dust thermal emis- 
sion, AGB stars are bright at infrared wavelengths and 
can be detected even towards highly obscured regions. 
Furthermore, the maser emission from their envelopes is 
strong enough to be detected throughout the Galaxy, and 
radio spectroscopic observations can provide th e stellar 
line-of-sight velocities to within a few kms~^ (e.g. lHabind 
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ously tar geted to search for 86 GHz SiO maser emission 
l)PaDer]|) . 

Here (Paper II) we present the available near- and mid- 
infrared photometry of the ta rgeted sources ( "SiO ta rgets" 
hereafter). In another paper ijMessineo et al.ir2003al here- 
after Paper III) we dea l with extinction corre ction and 
finally in the last paper l)Messineo et al.ll2003hl hereafter 
Paper IV) we compute and analyse the luminosities of 
the SiO targets. 

The spatial distribution of the 441 targets is shown in 
Fig.Q] These SiO targets are divided into two subsamples: 
253 sources were selected from the ISOGAL catalogue, 
("the ISOGAL sample"), and 188 sources from the MSX 
catalogue, ("the MSX sample"). The ISOGAL and MSX 
samples are examined to test whether both are drawn 
from the same "parent population" . Brightness variabil- 
ity is studied by a comparison of the DENIS and 2MASS 
photometry. 

The structure of the paper is as follows: in Sect. El 
we identify our SiO targets in various infrared catalogues 
and collect their magnitudes finding for many stars up to 
fourteen different measurements. In Sect. |31 we compare 
the statistical differences between our ISOGAL and MSX 
samples. In Sect. 0] we summarise additional information 
found with SIMBAD, e.g. variability and other types of 
masers, and we derive the probability of association be- 
tween the radio maser and the infrared counterpart. The 
brightness variability of the stars is discussed in Sect. [3 
In Sects. El and [71 we analyse the mid-infrared colours of 
the stars and compare them with those of OH/IR stars. 
The main conclusions are summarised in Sect.|Hl 

The individual source numbers (e.g #99) are taken 
from Table 2 (86 GH z SiO m aser detections) and Table 
3 (non-detections) in jPaoer j unless otherwise indicated. 
The SiO maser emission in this paper generally refers to 
the 86 GHz (w = l,J = 2^1) SiO maser only and not 
to the 43 GHz (J = 1 — 0) SiO masers. Velocities in this 
paper refer to line-of-sight velocities with respect to the 
Local Standard of Rest. 

2. Identification of the SiO targets in various 
infrared catalogues 

We cross identified all SiO targets, whether taken from 
the ISOGAL or from the MSX database, with all infrared 
catalogues available to us, the DENIS, 2MASS, ISOGAL, 
MSX and IRAS survey catalogues. The results are sum- 
marised in Table ^ In the following, we briefly recall 
the criteria used for the selection of the ISOGAL and 
the MSX samples and describe the modality of the cross- 
correlations. More details on the selection criteria can be 
found in .Paper L . 

2.1. ISOGAL data 

ISOGAL llOmont et al ll20n3^ is a 7 a,nd 1 5 ^m survey 
taken with ISOGAM l|Gesarskv et alJ Il996 f^ aboard the 
Infrared Space Observatory (ISO) satellite ()Kessler et al.i 



Il996l) . The 16 deg^ survey consists of selected sub- fields 
along the Galactic plane, mostly concentrating on the 
Galactic Centre. With a sensitivity approaching 10 mJy 
(two orders of magnitude deeper than IRAS and one order 
of magnitude deeper than MSX) and a resolution of 3-6" 
ISOGAL has detected over 100,000 objects. The ISOGAM 
data have been correlated with the DENIS near-infrared 
/, J and Ks band data to produce the five band ISOGAL- 
DENIS point source catalogue (ID-PSC) l|Schuller et alJ 
l2003il . The ID-PSC reports the photometric data in mag- 
nitudes. The astrometric accuracy of the ID-PSC is deter- 
mined by the present DENIS astrometric accuracy better 
than 0.5". 

Sources were selected from a preliminary version of the 
ISOGAL catalogue by their extinction-corrected 15 //m 
magnitude, [15]o, and their (i^so — [15]o) and ([7]o — [15]o) 
colours approximately corrected for extinction (Sect. 7.3.2 
jPaperll) . The brightest 15 /xm sources, [15]o < 1-0, and 
those with ([7]o-[15]o) < 0.7 and with (ifso-[15]o) < 1.95 
were excluded since they are likely to be foreground stars 
or non-AGB stars or AGB stars with very small mass-loss. 
Further, sources with [15]o > 3.4 were excluded since they 
are likely to show SiO maser emission fainter than our 
detection limit of 0.2 Jy. Sources with ([7]o — [15]o) > 2.3 
were excluded since they are likely to be compact HII 
regions or other young stellar objects or planetary nebulae. 
Those with (i^so — [15]o) > 4.85 were excluded because 
they are likely to be OH/IR stars with a high mass-loss 
rate or young stellar objects. Moreover, known OH/IR 
stars were discarded as the kinematic data are already 
known. 



2.2. MSX data 

The Midcourse Space Experiment (MSX) is a survey at 
five mid-IR bands ranging from 4.3 /im [Bl band] to 21.4 
fj,m [E band], with a sensitivity of 0.1 Jy in A band (8.2 8 
lim) and a spatial resolution of 18.3" I Price et'ani200lj) . 
The survey covers the Gala ctic plane to ±5 ° latitude. 
Version 1.2 of the MSX-PSC l|Egan et alJll999l) Hsts more 
than 300,000 point sources with an rms astrometric accu- 
racy of ~ 2". The MSX catalogue gives the source ffux 
density, F, in Jy. Magnitudes are obtained adopting the 
following zero-points: 58.49 Jy in A (8.26 /j,m) band, 26.51 
Jy in C (12.12 ^m) band, 18.29 Jy in D ( 14.65 um) b and 
and 8.80 Jy in E (21.41 ^m) band l|Egan et aljll999l) . 

For the MSX source selection we used flux densities 
in the A and D bands which have wavelength ranges 
roughly similar to the ISOGAL 7 and 15 fim bands. We 
selected those non-confused, good-quality sources in A 
and D band (flag > 3), which show variability in the 
A band. We avoided the reddest stars, Ed /Fa > 2.3 
{A — D > 2.2 mag). We also avoided the bluest and 
most luminous stars with Ed /Fa < 0.6 {A — D < 0.75 
mag) and Ed > 6 Jy (D < 1.2 mag) since they are 
likel y to be fore ground stars or supergiants 

(Schuller''20M 

■Schultheis et al...2003.) . Furthermore, following the classi- 
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Fig. 1. Location of the 441 SiO targets in Galactic coordinates. The 379 MSX counterparts are shown as open circles, 
the 267 ISOGAL counterparts as crosses. Overlap of MSX and ISOGAL sources resemple filled symbols. Four points 
fall outside the figure. 

Table 1. Number of counterparts of our SiO targets 
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fication of iKwok et al ] lll997l) of IRAS sources with low- 
resolution spectra, we used the C to E band ratio to dis- 
card very red {Fe/Fc > 1.4, C — E > 1.55 mag) sources, 
which are likely to be young stellar objects or OH/IR stars 
with thick envelopes. Known OH/IR stars were excluded. 

2.3. ISOGAL-MSX cross-identifications 

More than half (253) of our sample of SiO targets were 
selected from the ISOGAL survey, indeed from a prelim- 
inary version of the ID-PSC (Paper I). They are all in an 
area covered by the MSX survey. For each of our 253 ID- 
PSC selected SiO targets we searched for the closest as- 
sociated source in the MSX-PSC within 15" from the SiO 
target position. We found counterparts for 191 of them, 
190 of which were unique. The distribution of the angu- 
lar separations is shown in Fig. 12 the mean and median 
angular separation are 3.3" and 3.0" with a standard de- 
viation of 2.1", the maximum separation is 10.4". To es- 
timate the likely number of chance associations, we re- 



peated the cross-correlation after shifting all ID-PSC po- 
sitions between 30" and 100" finding that the probability 
of spurious associations within 11" is less than 1%. The 
62 ISOGAL sources without a counterpart in the MSX- 
PSC are mostly (55) concentrated in the central 3 degrees 
from the Galactic Centre. Their 7 /im ISOGAL fluxes 
range from 0.4 Jy to 1.5 Jy, whereas for the 191 associated 
sources, the flux is on average 1.6 Jy with a standard devi- 
ation of 1.3 Jy. Thus, the ISOGAL/non-MSX sources are 
likely to have been excluded from the MSX-PSC because 
of confusion due to the high stellar density in the inner 
Galaxy. The fraction of our ISOGAL stars without MSX 
identifications is ~11% for those with 7 /itm magnitude 
[7] < 4.5 (^^7 > 1.3 Jy). 

Out of the 188 SiO targets selected from the MSX 
catalogue, only 14 are located in an ISOGAL field, and 
we found ID-PSC identifications for all of those within 5" 
from the MSX positions. For the remaining 174 sources 
we searched the DENIS PSC. 
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Fig. 2. Associations between our SiO ISOGAL targets 
and the MSX catalogues. The distribution of the angular 
separations (") between the ISOGAL and the MSX posi- 
tions (continuum line). For comparison, the dashed line is 
the total distribution of the chance associations, found by 
performing five shifting experiments. 



2.4.1. ISOGAL sample and DENIS identifications 

For the SiO targets selected from the ISOGAL survey, 
a proper DENIS iden tification and photometry had a l- 
ready been provided l|Schuller et al.lbnO,'^ ISimonI 1200.'^ . 
For each ISOGAL field, a sample of DENIS point sources 
from the same area was extracted. DENIS coordinates 
were retained as the astrometric reference system and 
a polynomial distortion correction was applied to the 
ISOGAL coordinates in order to match as best as pos- 
sible the ISOGAL and DENIS reference coordinates. In 
order to reduce the number of spurious associations, only 
DENIS sources with a Ks detection were selected and a 
JTr upper limit was imposed according to field density 
ijSchuller et al.ll2003 l). Among our 253 ISOGAL targets all 
have a DENIS Ks identification, while only 86% and 21% 
have J and / detections, respectively. 



Fig. 3. Difference between the MSX band D magni- 
tude, D, and the ISOCAM 15 fiui magnitudes, [15], versus 
the D. Mag nitude zero point in D band was taken from 
lEgan et aTTIlQQfll . The continuous line is I? - [15] 
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As the ISOCAM 15 ^m and MSX D filters are simi- 
lar, we compared the ISOCAM 15 /im magnitudes, [15], 
and the MSX D band magnitudes, D, of the 154 sources 
detected at 15 /zm in both surveys and found good agree- 
ment (Fig. 131; the average difference D — [15] is —0.04 
magnitude and the standard deviation is 0.3 mag, result- 
ing from the combination (~ 0.15 mag) of the photometric 
errors of both catalogues and from the possible intrinsic 
source variability. 

2.4. DENIS data 

DENIS is a simultaneous / (0.8 fim), J (1.25 /im), and 
Kg (2.15 /im) band survey of the southern hemisphere us- 
ing th e ESO 1-meter telescope at La Silla ijEpchtein et alJ 
Il994j) . The 3ct detection limit in the respective bands is 
0.05 mJy (19 mag), 0.5 mJy (16 mag) and 2.5 mJy (13.5 
mag), and the saturation limit is at magnitudes 10, 7.5 
and 6, respectively. The absolute DENIS astrometry is 
fixed using the USNO-A2.0 catalogue. The current abso- 
lute astrometric accurac y of DEN IS is better than 0.5" 
(rms) (see Sect. 4.2 in iSchuller et al...2003i) . 



2.4.2. MSX sample and DENIS identifications 

We searched for possible DENIS counterparts to the 188 
SiO targets which we had selected from the MSX cata- 
logue. For all sources from the MSX sample we searched 
the provisional bulge DENIS PSC dSchnller et a].ll20fl3t 
Isimon 200 3) for the nearest Kg band object. As late-type 
stars are intrinsically red, only Ko, band counterparts were 
examined. As saturated sources are usually given in the 
catalogue with null magnitude, DENIS p hotometry from 
other observations ((Epchtein et al.iri994l) and the Ks im- 
ages were used as additional checks for saturated counter- 
parts to the MSX SiO targets. Fortunately, such DENIS 
data was available for all MSX sources that we observed; 
however, a large number of bright sources saturated the 
detector in the Ks band. While all sources were detected 
in the Ks band, the fraction of detections in the J and / 
band was 79% and 22%, respectively. 

The distribution of angular separations of the Ks band 
counterparts identified is shown in Fig. 2J The mean and 
median separations are 3.0" and 2.7" with a standard de- 
viation of 1.8", respectively, which is consistent with the 
expected scatter due to positional uncertainties for the 
MSX and DENIS catalogues. 

To find the distribution of random chance associations 
we again searched the nearest neighbours after shifting 
the coordinates of the SiO targets by 30". The result- 
ing distribution of separations and magnitudes is also 
shown in Fig. 01 The real associations are on average 
much brighter and closer than the chance associations. To 
compute the expected number of incorrect identifications 
we divided the "chance" distribution in magnitude bins, 
(m,TO -|- Am), and in each bin computed the normalised 
cumulative chance distribution of the separations, f{m, r), 
i.e., the fraction of all chance associations with those mag- 
nitudes and within a radius r. We then computed the sum 
I]/(m,r) over all "real" associations, each being charac- 
terised by m and r. Thus, f{m, r) gives the probability for 
an association to be spurious, and the sum over all sources 
yields the total expected number of spurious identifica- 
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Fig. 4. Associations of our SiO MSX targets in the 
DENIS catalogue. Top panel: distribution of Ks mag- 
nitudes for the associations (continuum line) and for the 
chance associations (dashed line). Lower panel: The con- 
tinuous line shows the distribution of angular separations. 
The dashed lines show the normalised cumulative distribu- 
tions, f(r,m), of chance associations found by misalign- 
ing the catalogue and MSX SiO targets, where we plot 
separately, starting from the bottom, the distributions of 
sources with i^s < 9, 9 < i^s < H, and Ks > 11 mag, 
respectively. 



tions. When we consider only the brightest identifications 
for the 154 sources with a i^s < 9 mag counterpart, only 
two spurious identifications are expected. For the 34 pos- 
sible identifications with Kg > 9 mag, we would expect 
three to be spurious. In total, we expect that about five 
of our identifications may be incorrect. 

Since for the fainter possible counterparts the chance 
of a false identification is higher, we looked for brighter 
possible identifications somewhat further away, which due 
to the lower surface density of brighter sources might have 
a higher probability to be the actual counterpart. In a 
few cases (#231, #243, #367, #405, #406, and #424) 
we found brighter sources somewhat further away but 
with a lower value of f(m,r) than the closest identifi- 
cation. This suggested that these brighter sources were 
more likely to be the correct counterparts, which we 
therefore retained. For these MSX targets with dubious 
near-infrared association, we additionally examined other 
MSX sources in their surrounding and checked for possi- 
ble astrometric shifts between the DENIS and MSX co- 
ordinates which could uniquely identify the correct near- 
infrared counterpart. However, because of the low MSX 
source density 1 source per 4'x 4'), only few asso- 
ciations could be confirmed in this way. Seven sources 



could only be associated with a. Kg > 11 mag counter- 
part (#162,#273,#363,#403,#417,#423,#443) within 
15" (beam size of the IRAM telescope). 

2.5. 2MASS data 

The Two Micron All Sky Survey (2MASS) has surveyed 
the entire sky at near-infrared wavelengths from the 
Whipple Observatory on Mt. Hopkins, AZ, and the Cerro 
Tololo InterAmerican Observatory (CTIO). The cameras 
observed simultaneously in the J (1.25 /xm), H (1.65 fim) 
and Kg (2.17 fim) bands with a 10 cr sensitivity limit in 
uncrowde d fields of 15.8, 15.1 and 14.3 mag, respectively 
(jBcichnia n et al Jl998tlCutri et al.l200^ . The absolute as- 
trometry has typical uncertainty of 0.1" (rms) and is based 
on the Tycho-2 and UCACrlO catalogue (a new version of 
the USNO's ACT catalogue). The on-line 2MASS PSC (aU 
sky data release) is accessible at the Infrared Processing 
and Analysis Center (IPAC). 2MASS data was available 
for all our positions. We retrieved 2MASS sources for our 
441 SiO target positions and we searched for the closest 
positional match. Since the astrometric accuracy is bet- 
ter than 0.5" for the ISOGAL targets but only -2" for 
the MSX targets, a different distribution of angular sep- 
arations is expected for the MSX-2MASS and ISOGAL- 
2MASS associations. 
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Fig. 5. 2MASS Ks image 230" x 230" centred at the po- 
sition of #164. The black circle delimits the main beam 
(FWHM) of the IRAM telescope used for the SiO maser 
search. A small gray circle marks the 2MASS association 
and finally crosses mark sources detected by ISOGAL. In 
agreement with the discussion in the text, all the ISOGAL 
and 2MASS positions are in excellent agreement. Note 
that the targeted source is the only mid-infrared source 
falling inside the IRAM telescope beam. 
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To avoid misidentification, because of the high source 
density of the 2MASS survey and because the DENIS 
counterparts are mostly brighter than 11 mag in Kg, we 
limited the 2MASS Kg magnitude to Ks < 11. However, 
seven sources in our sample, #162, #273, #363, #403, 
#413, #421 and #423, could be associated only with 
2MASS sources fainter than ifg = 11 (those sources are 
also faint in DENIS) within 15", which is the beam size 
of the IRAM telescope. Positional associations were con- 
firmed via overplotting the 2MASS counterpart image 
with both the SiO targets and the DENIS sources. Finding 
charts were obtained for all the stars with 2MASS im- 
ages, an example of which is given in Fig. (SJ We found 
439 2MASS counterparts and missed only two. In fact, af- 
ter image inspections, the potential 2MASS counterparts 
for two sources, #224 and #298, were eliminated as their 
positions on the 2MASS images were marked by artifacts. 

2.5.1. ISOGAL sample and 2MASS identifications 

The mean, median sep arations between the ISOGAL SiO 
targets (positions as in IPaoer^ . without rounding RA to 
one tenth of second) and the 2MASS associations are 0.7" 
and 0.3" with a standard deviation of 1.0", (see Fig.^l. 

There is a non-Gaussian tail at large separations. We 
have individually checked all the sources with separation 
larger than 3" and note that they all have Kg < 6. We 
attribute the large separation to saturation of the DENIS 
detector. Saturated pixels are an obstacle to the correct 
determination of the source centroid and this affects the 
astrometry of saturated stars. Furthermore, most of these 
bright sources do not have any / associations and therefore 
the J / Kg astrometry is kept. The mean, median and stan- 
dard deviation of the separations between the ISOGAL 
SiO targets and the 2MASS associations with Ks> 6.5 
mag are 0.4", 0.3" and 0.4"; while for 2MASS associations 
with Ks< 6.5 mag they are 1.9", 2.1" and 1.5". 
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Fig. 6. Distribution of the an gular sep arations between 
the SiO targets (positions as irL lPaperl) and the 2MASS 
associations. The continuum line shows the distribution 
of the ISOGAL SiO targets and the corresponding y-axis 
is on the left side. The dashed line shows the distribution 
of the MSX SiO targets and the relative y-axis is on the 
right side. 



2.5.2. MSX sample and 2MASS identifications 

We found 187 2MASS counterparts of MSX SiO targets. 
The mean, median and standard deviations of the sep- 
arations between the MSX SiO targets (using the MSX 
coordinates) and the 2MASS associations are 2.9", 2.6" 
and 1.8", the same as the separations between the MSX 
SiO targets and the DENIS associations (see Sect. I2.4.'2|l . 
This is again consistent with the expected scatter due to 
positional uncertainties of sources in the MSX and 2MASS 
catalogues. The mean, median and standard deviation of 
the separations between the MSX SiO targets (using the 
DENIS coordinates) and the 2MASS associations are 0.6", 
0.4" and 0.7", respectively. The mean, median and stan- 
dard deviation of the separations in right ascension RA 
are 0.0", 0.0" and 0.5"; while the mean, median and stan- 
dard deviation of the separations in declination DEC are 
-0.2", -0.1" and 0.6". 

2.6. SiO targets and IRAS identifications 

We checked for counterparts to our SiO targets in the 
Infrared Astronomical Satellite Point Source Catalogue 
(IRAS PSC). Toward the galactic plane the IRAS sur- 
vey is strongly limited by confusion at the low resolution 
of the IRAS instruments (0.5' at 12/im), especially in the 
longer wavelength bands. Since much of the past work on 
maser surveys of late- type stars is based on the IRAS data, 
it may be useful to have an according comparison. 

We selected only the 165 IRAS associations within 3.5c7 
error ellipse of the IRAS PSC, to reduce the chances of 
spurious associations to ^ 2% (3 sources). A compari- 
son of 12/im fluxes of the prospective counterparts with 
the ISOGAL and MSX fluxes (Fig.|Hll shows a good agree- 
ment, confirming that the IRAS identifications are proper. 
About 35% of our SiO targets have counterparts in the 
IRAS PSC: 65% of the MSX sample and 17% of the 
ISOGAL sample. Of those, 96% are detected at 12/im, 
87% at 25/im, 6% at 60/im and 4% at 100/im, and 56% 
are reported in the IRAS catalogue as variables (flag> 80). 

2. 7. SiO targets and visual identifications 

Within a radius of 5" to our SiO targets we searched for 
possible visual counterparts in the Tycho 2 and USNO- 
A2.0 catalogues. We found possible associations for 85 
SiO targets, of which 27 are validated by corresponding 
/ band counterparts. Their R magnitudes range from 10.9 
to 17.9 and B magnitudes from 13.4 to 20.6. Considering 
that the pulsation amplitude of Mira stars increases at 
short er wayeleng hts and can be up to 8 mag in the vi- 
sual |2ma3^^J, all those visual stars are possible coun- 
terparts of our SiO targets. Most of them are located 
at latitude |5| < 0.8°, but the extinction value inferred 
by their colours are mu ch smaller than the median of 
their surrounding stars ijPaper II J) . Therefore they are 
likely to be foreground stars. Two SiO targets, #7 with 
{l,b,vel) = (-1.50°, 0.95°, 18.9 kms'^) and #139 with 
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Fig. 7. Upper panel: Separation of the possible IRAS 
counterparts, expressed in IRAS sigma units, against an- 
gular separation. Crosses indicate IRAS sources with pos- 
sible MSX counterparts closer than those associated with 
the SiO targets; these are likely to be unrelated to the 
SiO targets. The dotted horizontal line is the upper limit 
that is selected. Lower panel: The continuum line shows 
the distribution of the separations (in sigma units) of the 
possible IRAS counterparts. The dashed line shows the 
distribution of the chance associations obtained shifting 
the source coordinates by 250". The dotted vertical line is 
our chosen upper limit of 3.5 a. 



{l,b,vel) = (0.31°, -2.18°, 205.7 kms^^), have extinction 
value inferred by their colours A^— 5.5 and 3 mag, re- 
spectively, consistent with the median extinction of their 
surrounding stars (PaDcrlll). Furthermore, the velocity 
of #139 is inconsistent with being a foreground star. 
Therefore we conclude that they are likely located in the 
Galactic bulge, in regions of low interstellar extinction. In 
fact, #139 is located in the optical window WO. 2 — 2 1 
at (Z,&) = (0.2°, -2.15°) llDutra et al.ll2002l:IStaneklll998l) . 
and #7 has an exti nction value typic al of bulge ISOGAL 
fields with 5 - +1° l)Oiha et al.ll2003|) . 

2.8. The table 

Tabic 121 lists the infrared photometry of the SiO targets. 
The columns of Table |21 are as follows: an identification 
number (ID), the same as in Table 2 and 3 of IPaoer"^ 
followed by the Right Ascension (RA), and Declination 
(DEC), in J2000 of the 2MASS counterpart; the DENIS 
(/, J, Ks) magnitudes, the ISOGAL ([7], [15]) magnitudes, 
the 2MASS (J, H, Ks) magnitudes; the angular separa- 
tion (dis), between the adopted near-infrared position and 



X 

— J^^^.^v 




1 ' ' ■ ' IX 

X 

* y ' 
w ^x X 






1 < < . < 



[121 

Fig. 8. Difference between the MSX magnitude in the 
C band, and the IRAS 12/im magnitudes, [12], (filled cir- 
cles) versus the [12] . The 18 open circles indicate sources 
without C measurement, for which the 15yLtm ISOGAL 
measurement is plotted. Crosses are likely spurious MSX- 
IRAS associations (angular separation larger than 3.5ct). 
The plotted IRAS magnitudes are -2.5 log(i^[Jy]/28.3). 
The continuous line is C— [12] = 0.00 mag and the dashed 
line is C — [12] = 0.05 mag, which is the mean C — [12] 
colour. 

the MSX position, the MSX {A, C, D, E) magnitudes, the 
IRAS 12 and 25 /im magnitudes and finally a variability 
fiag (war), defined as described in Sect. [3 An additional 
column is used for comments on individual stars. 



3. A comparison between the ISOGAL and MSX 
samples 

A large variety of names exists to indicate oxygen-rich 
AGB stars characterized by different pulsation properties 
and/or mass-loss rate: semi-regular (SR) stars and Mira 
stars (Ha in emission, visual pulsation amplitude larger 
than 2.5 mag), large amplitude variables (LAV), long pe- 
riod variable (LPV) stars (when their periods are longer 
than 100 days), and OH/IR stars (with 1612MHz OH 
maser emission). In the IRAS colour-colour diagram the 
oxygen-rich AGB stars are distributed on a well-defined 
sequence of increasing shell opacity and stellar mass-loss 
rate (e.g. Habing 1996), which goes from SRs and Miras 
with the bluest late-type colours and the 9.7 /im silicate 
feature in emission, to the coldest OH/IR stars with the 
reddest colours and the 9.7 /im silicate feature in ab- 
sorption. The sequence of increasing shell opacity corre- 
spo nds also to an increasi ng JTs- [15] or JTr- [12] colour 
fe.g.lOiha et al...2003; .Olivier et alJl200lHWhitelock et alJ 



11994 Paper~ 



SiO maser emission is generated in the envelopes 
of mass-losing AGB stars, close to their stellar photo- 
spheres and it occurs more frequently towards oxygen-rich 
Mira stars than towards other AGB stars (includ ing SR 
and OH/IR stars) (Buiarrabal 1994:lNvman et aL 1993) . 
Therefore, for our 86 GHz SiO maser survey we selected 
the brightest sources at 15 /xm with colours of Mira-like 
stars (fPaogrl ). The ISOGAL and MSX samples were se- 



Table 2. Infrared counterparts of the SiO targets.* The identification numbers are the same as in Table 2 and Table 3 of lPaper 
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* The full table is available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb.u — strasbg.fr/cgi — bin/qcat?J/A -|- A/ (vol) /(page). 
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[7]-[15] 

Fig. 9. Mid-infrared colour-magnitude diagrams. Left-hand panel: all data from ISOGAL. Right-hand panel: 

all data from MSX. Filled circles indicate targets from the ISOGAL sample and crosses indicate targets from the 
MSX sample. The two samples have similar mid-infrared colours, due to selection criteria. The numbers of MSX and 
ISOGAL targets between the two panels vary as explained in Sect. 12.31 



lected to have similar 7 and 15 /im colours, as shown 
in Fig. 1^ For the ISOGAL sample a range of intrinsic 
(i^'s — [15])— colour was selected (see Sect. 12. If) , but for the 
MSX sample no near-infrared counterparts were available 
at the time of the observations, and thus no Kg magni- 
tudes. As an alternative to the (Tis — [15]) criterion, for 
the MSX sample we imposed an upper limit to the ra- 
tio of the fluxes in the E (21/im) and C (12/^m) bands 
{C — E < 1.55 mag). Both criteria were defined in order 
to avoid non-variable objects and thick envelope objects, 
but these criteria are not equivalent. The emission in the 
Ks band is dominated by the stellar emission attenuated 
by the circumstellar absorption, while circumstellar dust 
emission contributes strongly to the mid-infrared radia- 
tion. 

AGB stars with thin envelopes have C — E <^ 1.5 
mag, while sources with C — E >~ 1.5 mag are AGB 
stars with thick envelopes , post-A GB stars and young 
stars l|Lumsden et al.ll2n02t FSevenster .20Q2) . Our sample 
of SiO targets includes very few sources with C — E > 1.5. 

The 7is~[15] and i^g — [12] are good indicators of mass- 
loss rate for AGB star with shells at few hundred Kelvin. 
All Miras have redder i^s — [12] than non-Mira stars from 
1.8 mag up to 6 or ev en 14 mag (e.g. lOlivier et al.ll2n0lt 
IWhitelock et alJll994l) . Thick-envelope OH /IR stars have 
typically Ks-[15] > 4 llOrtiz et a] .1120021) . The ISOGAL 
SiO target sample is characterised by a smaller range of 
(iis — [15]) or (Ks—D) values than the MSX sample, as 
shown in Fig. ^1 suggesting that the MSX sample in- 
cludes a tail of sources with optically thick envelopes. 
How ever, this will be ver ified after correction for extinc- 
tion ijPaDerlllL IPaperlV^ . 



Figure [TUl shows a strong correlation, between Kg and 
{Ks — [15]) or (Ks—D), which is due to the way we selected 
our original sample. In fact, due to the general correlation 
of SiO maser emission and infrared luminosity (Bujarrabal 
et al. 1987), we selected only sources with [15] < 3.4 in 
order to be able to detect the expected SiO line. The [15] 
and D magnitudes range from 3.4 to 1.0. This narrow 
range generates the correlation seen in Fig. 1101 

Figure 10 also shows that the two samples, ISOGAL 
and MSX, overlap largely, although there are some minor 
systematic differences. There is a vertical shift between 
the ISOGAL and MSX sequence, with the MSX sources 
brighter for a given colour. This effect is due to the dif- 
ferent sensitivity between the MSX D band and ISOGAL 
surveys. MSX targets have on average a brighter D (or 
[15]) magnitude than the ISOGAL targets (see Fig. 
This fact translates in differences in distance between the 
two samples and suggest that the MSX stars are closer on 
average; the two samples are also distributed differently 
in longitude (see Fig. 

On the basis of this comparison in the following we 
combine the results obtained with the ISOGAL sample 
with those obtained with the MSX sample, taking into 
account that there are some minor differences between the 
two samples. 



4. Some other remarks on the SiO targets 

4.1. SIM BAD search 

A SIMBAD search revealed some extra information about 
our SiO targets. 
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[15] = 3.4\.. 



[15]=1 



[K3]-[15] 

Fig. 10. 2MASS Ks versus {Ks-[15]) or {Ks-D). Filled 
circles indicate objects from the ISOGAL sample, for 
which we plot (ii's — [15]); and crosses indicate objects from 
the MSX sample, for which we plot (Ks—D). The two 
dashed lines correspond to [15] = 3.4 and 1.0 mag. 



Sources #7, #286 and #303 are included in the cat- 
alogue of late-type s tars in the inner Galactic region by 
iRaharto et"al] l|l984l) as spectral types M6, M7 and M6.5, 
respectively. 

#153 is a well known Mira star, TLE 53 , with a pe- 
riod o f 480 d, located in Baade's window (e.g. lGlass et alJ 
19951). Our samp le also includes 15 LPVs found by 
Glass et al.l ()200ll) within 0.3° from the Galactic cen- 
ter and 19 candidate variable stars from the list of 



difference between the SiO heliocentric velocity Vhei — 
+4.4 ± 1.0 kms^^ and th e optical hel iocentric velocity 
Vhei = -14.5 ± 1.0 kms-i (|Torrealll999l de La Reza priv. 
communication). Since the SiO maser velocity is usually 
coincident wi th the stellar velocity within few kms~^ (e.g. 
lHabindll993) . we suggest that the SiO emitter is not as- 
sociated with the G8II star, IRAS 17590-2412/PDS 482, 
which however is the only mid-infrared source within the 
IRAM beam and the association between the ISOGAL 
and DENIS source is of excellent quality (flag=5). 

#189, IRAS 18059-2554, is given by 

iLvnch fc Rossanol l)l990|) as a possible member of 
the globular cluster NGC6553. The stellar line-of-sight 
velocity, obtained through the SiO maser line, is 161.2 
kms^^. Because the cluster mea n line-of-sight veloc ity is 
7 kms"^ with a = 
conclude that IRAS 
cluster. 

Twenty-eight of our 86 GHz SiO targets were previ- 
ously observed for 43 G Hz SiO maser emission. These are 
discussed in Sect. 4.5 o f lPaperj . For the strongest 86 GHz 
maser sources within 2.2° of the Galactic Centre from 
PaperJ we recently used the Very Large Array (VLA) to 
observe the two 43 GHz SiO maser line s (v=:0 and v=l) 
simultaneously l|Siouwerman et alllioosfl . 

We excluded from our SiO maser s ur vey the 
OH/I R stars detected b y ISeyenster et al.l lll997alhL 
1200 ih ISiouwerman et all l|l998l) and iLindavist et alJ 



14 kms"^ (jCoelho et al.ll200l[l . we 

18059—2554 is not a member of the 



1119921) . However, 
and the limited 
surveys, we still 
with a SIMBAD 



due to intrinsic source variability 
sensitivity of the Sevenster et al. 
included 4 OH/IR stars, as found 
search: #181, #226 and #257, aU 



ISchtiltheis et alJ ||200G ') (listed in Table 2 and 3 of lPaperJ) . 

A few sources are given in the literature as possible red 
supergiants or extremely luminous AGB stars: #25, #32, 
#92 an d #295 c orrespond to sources #6, #8, #31 and #5 
of Nagata et al.l (il993), respectiv e ly; #3 56 is classified as 
bulge M supergiant by iRahartol l|l99l[) and IStephensonI 
lll992l) also lists it among distant luminous early type 
stars. 

#127 (IRAS 17500-2512), #178 (IRAS 18040-2028), 
#188 (IRAS 18060-1857), #252 (IRAS 18285-1033), 
#265 (IRAS 18367-0507) and # 434 (IRAS 18415-0355) 
are listed by ijKwok et al.lll997|) among sources detected 
with the IRAS Low Resolution Spectrometer, in the range 
8-23 /im and with a resolution, A/AA, ~ 20 — 40. Four 
spectra are noisy or incomplete, while the spectra of #178 
(IRAS 18040-2028), #434 (IRAS 18415-0355) are clas- 
sified as featureless; they are probably evolved stars with 
negligible amounts of circumstellar dust. For those two 
objects, at longitudes 9.7 and 28.6° respectively, we also 
comp ute a mod erate mass- loss rate of 5 - 7 x lO"'^ Mq 
yr~i l|Paper lit) . 

#164, IRAS 17590-2412 is cl assified as a Li K gi- 
ant star bv lde La Reza et a,1.l \^99i(\ . There is a significant ^ 



with detected SiO maser emission, coincide in posi- 
tion and velocity with 0H9. 84+0.01 , OH17.43-0.0 8 
and OH25. 05+0.28, respectively (Blom maert et alJll994|) : 
#409 (IRAS 18142-1600), not detected in our SiO sur- 
vey, corresponds to OH #280 (OH14.805+0.150) listed by 
Ire Lintel Hekkert et al.. (19891 



4.2. Is the targeted star the actual SiO emitter? 

Using the IRAM-30m telescope, we detected 86 GHz 
SiO maser emission toward 268 of our targeted posi- 
tions, where three positions showed a double detection 
(#21-#22, #64-#65 and #78-#79). ^ To confirm that 
our targets are the actual SiO emitters, we need to ex- 
amine all other possible objects falling inside the 29" 
(FWHM) IRAM beam. Due to the correlation between 
intensity of the maser line and mid-infrared brightness 
(|Buiarrabal et al]ll987j) . we can restrict the analysis to 



^ Throughout the paper for statistic computations only one 
infrared source, the targeted, is considered for these three posi- 
tions. A second mid-infrared source is found within the IRAM 
beam at the position of the SiO sources #78 and #79 (sepa- 
ration ^ 12"), and at the positi on of th e SiO sources #21 and 
#22 (separation of ^ 15") (see lPaperll) and is given in Table 
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only mid-infrared objects. Considering the 379 SiO tar- 
gets with an MSX identification, we found in all but one 
case only one target inside the 86 GHz beam. Considering 
the 267 SiO targets with an ISOGAL identification, 89% 
are unique ISOGAL objects inside the main beam. SiO 
maser emission was detected towards 19 of our 29 posi- 
tions with more than one ISOGAL source in the beam 
and in all cases the targeted star is the brightest 15 /zm 
source. 
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Fig. 11. Upper Panel: Difference between 2MASS and 
DENIS J magnitudes versus the DENIS J magnitude. 
Filled squares represent the ISOGAL SiO targets. For 
comparison small dots show ID-PSC and 2MASS asso- 
ciations obtained in several ISOGAL fields, which have a 
distribution of their J variations consistent with a gaus- 
sian distribution centred at zero. No correction for offset 
in the photometric zeropoint was applied. Open circles 
represent the MSX targets, which have indications of vari- 
ability at mid-infrared wavelengths. The error bars shown 
within the box arc ±1(7 mag. The vertical line indicates the 
DENIS saturation limit. Lower Panel: Comparison of 
the 2MASS and DENIS Kg magnitudes versus the DENIS 
Kg magnitude. Symbols are as in the top panel. 



5. Variability 

The 86 GHz SiO maser line intensity is observed to be 
stron ger in 0-rich Mira stars than in other types of AGB 
stars (jBuiarraball 11994 iNvman et alJll99i) . To increase 
our chance of detecting the SiO maser emission, we tried 
to select strongly variable sources ( Paper I, and references 
therein). Therefore, our selected MSX targets all have an 
indication of variability in the A band. Unfortunately, 
the ISOGAL database does not provide direct informa- 



tion on variability. However, from their position in the 
ISOGAL/DENIS {Kso - [15]o vs. [15]o) diagram we were 
expecting that our ISOGAL SiO targets were mainly large 
amphtude variables (Sect. 4.2 of IPaperJ^ . For 62 ISOGAL 
targets out of the 191 which have an MSX counterpart, 
there is an indication of variability in the A band (or for 
94 in at least one of the A, C, D, and E bands). 

5.1. Variability information from DENIS and 2MASS 
data 

The comparison between the near-infrared photometry 
obtained during the course of the DENIS survey and that 
of 2MASS again confirms our hypothesis for variability in 
the ISOGAL sample. 

The J and Kg filters used by DENIS and 2MASS are 
similar and therefore the measurements obtained during 
the course of the DENIS and 2MASS surveys are directly 
comparable. For non variable sources, the differences be- 
tween the J magnitude of DENIS and 2MASS and the Kg 
mag nitude of DENIS and 2MASS are smaller than 0.15 
mag dPelmotte et aLlEfflilSchultheis G]a,SFJl2nmL and 



present work). 

The DENIS observations were performed between 1996 
and 2000, while the relevant 2MASS observations were 
performed between 1998 and 2000; AGB variables have 
periods from 50 to 1000 days, therefore the interval of 
time between the observations makes it possible to derive 
variability information. 

For each of the 61 ISOGAL fields containing our 
SiO targets we retrieved the corresponding 2MASS sub- 
catalogue, and cross-correlated the ID-PSC and the 
2MASS point source positions. For our ISOGAL sample 
the difference between the 2MASS and DENIS J and Ks 
magnitudes is shown as a function of the corresponding 
DENIS magnitude in Fig. ^ Figure ^ shows the differ- 
ence between the 2MASS and DENIS J magnitudes plot- 
ted against the time between the 2MASS and DENIS ob- 
servations. The distribution of the magnitude variations 
of the ISOGAL targets is different compared to that of 
random field objects. The Kolmogorov-Smirnov test gives 
a zero probability for the ISOGAL targets and field stars 
to be extracted from the same population. For 55% of our 
ISOGAL stars, the difference between both the 2MASS 
and the DENIS J and the 2MASS and the DENIS Ks 
magnitudes is larger than 3 times the dispersion measured 
in the corresponding field. Therefore, our sample contains 
mostly variable stars. 

Due to the simultaneity of the J and Kg measurements 
in both the DENIS and 2MASS surveys, a correlation is 
expected between the variation in the J magnitude (A J) 
and in the Kg magnitude (Aiig). As shown in Fig. 1 131 such 
a correlation exists. A linear least squares fit yields 

AJ = 1.57(±0.03) X AKs + 0.05 ± (0.01). 

The relative pulsation amplitude in the Kg band is ^ 60% 
of the relative amplitude in J band. 



12 



M. Messineo et al.: 86 GHz SiO maser survey of late-type stars in the Inner Galaxy 



Z 

w 

Q 



ITl 
< 




Fig. 12. The 

nitude versus 
observations. 



100 200 300 400 500 
Time difference [days] 

difference of the 2MASS and DENIS J mag- 
the time between the 2MASS and DENIS 



S 2 



SiO largel.s 




Field slrti.- 



-2 2 

Kj {2MASS) - Kg (DENIS) 



-2 2 

Kg (2MASS) - Kg (DENIS) 



Fig. 13. Difference between the 2MASS and DENIS 
J magnitudes versus the difference of the 2MASS and 
DENIS i^s magnitudes. Left panel: SiO maser targets. 
Upper limit measurements and measurements affected by 
saturation are excluded. The photometry of #42 is af- 
fected by blending with another source. The continuous 
line shows our best fit. Within the box the typical error 
for sources in our range of magnitude is drawn. Right 
panel: for comparison, field stars with J < 14 mag. 



For comparison, in Fig. we also show the rela- 
tion between the pulsation amplitudes in the J and K 
SAAO bands for two different samples of oxygen- rich stars 
in the solar neig hbourhood llOlivier et al."'200l') and the 
South Galactic Cap (SGC) (IWhitelock et al...l99l . The 
two samples have a different period distribution; most of 
the stars from the solar neighbourhood sample have peri- 
ods between 500 and 700 d, while most of the SGC stars 
have periods between 150 and 450 d. Overplotting our 
best-fit, we see that it aligns well with the distribution of 
the two samples of LPVs. 

A monitoring program of the near-infrared magnitudes 
of our SiO maser sample will provide pulsation periods 
and estimates of the source distances through the period- 
luminosity relation. 



5.2. Variabity flag 

All the available information on variability, from MSX 
data (variability in at least one of the MSX band) , IRAS 
data (variability flag > 50), DENIS- 2MASS data (AJ > 
3cr(J) and AK > 3a(K )) or f rom iGlass et all ll200lh 
from ISchultheis &: Glasd l)200l|) and from IClass et alJ 
l|l995l) , is summarised defining a fiag (var) which is 2 when 
variability is detected in at least one of the datasets, 1 
when variability is not detected but not all the datasets 
were available and when variability is not detected in 
any of the datasets. All the 188 MSX targets have the 
variability flag equal to 2 (due to selection). Among the 
253 ISOGAL targets 150 have the flag equal to 2, 81 equal 
to 1 and 22 equal to 0. 

There is no correlation between the variability indica- 
tion and the detection of SiO maser emission. 




Fig. 14. Pulsation amplitude in J band versus amplitude 
in K band for two samples of oxygen-r ich variable stars : 
one in the solar neighbourhood from Olivier et alJ l|200j) 
(star red points) and a.nothe r in the South Galactic Cap 
from lWhitelock et a l.* 1*1994^ (filled circles). Superimposed 
is our best fit from Fig. E| 



6. The distribution of the SiO targets in the IRAS 
two-colour diagram 

We cannot determine the distribution of our tar- 
get sources in the t wo-co lour IRAS diagram of 
Ivan der Veen fc Habind l)l988j) since for most of our 
sources we have only upper flux limits at 60 /im. Figure 
115! shows the distribution of the IRAS [12] — [25] colours 
and recalls the regions of the IRAS two-colour diagram 
that separate different classes of evolved stars with cir- 
cumstellar envelopes, from the bluer Mira to the thick en- 
velope OH/IR stars. The [12] — [25] colours of our selected 
sources range from —1.0 to 0.4, peaking at —0.2, corre- 
sponding mostly to regions Ilia and VIb. Region Ilia rep- 
resents sources with moderate dust emission, being pop- 
u lated m ostly by oxy gen-rich st ars with silicate emission 
(Kwok et al. 1 997t ivan der Veen fc Hab ing 1988). Region 
VIb, whith 60/im excess, contains a mixed population of 
early type stars with line emission and planetary nebula; 
however this region is scarcely polulated. Considering the 
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Fig. 15. Distribution of the IRAS [12] - [25] colours, 
— 2.5 log(i^i2/F25)- The continuous thin line indicate the 
distribution of our SiO targets; the thick line indicates 
the distribution of SiO targets with Kg < 6.5 mag. For 
comparison, the dashed thin line indi cates the d istribu- 
tion of the colours of OH/IR stars of ISevensteil l|2002l) . 
Upper limit flux densities are excluded. In the upper right 
corner, we recall the IRAS two-colour diagram. 



distribution of the IRAS good quality sources, only 5% of 
them are located in region VIb. We conclude that our tar- 
gets are mostly Mira stars with moderate mass-loss rate, in 
agreement with the selection criteria. For comparison, the 
distribution of the [12] — [25] colours of the OH/IR stars of 
ISevenster {.2QQ2) is also shown in Fig. El OH/IR stars are 
distributed over a larger and redder colour range, although 
partially overlapping with the colours of the SiO sample; 
they can have significantly (~ 1 mag) redder colours than 
the SiO targets. This can not be accounted for by interstel- 
lar reddening, which is only 0.1-0.2 mag for Ay ~ 20 — 30 
mag. The stars in our sample brighter than Kg < 6.5 mag 
are bluer than those with Kg > 6.5 mag. They are likely 
to be mostly foreground stars with thinner shell, but still 
IRAS detectable due to their proximity. 

#177 (IRAS 18039-2052) is the only source with 
[12] — [25] > 1.0, which resembles those colours of an 
embedded young stellar object. However, methanol and 
water maser emissi on was unsuccessfully s earched for to- 
war d this source (Mac Leod et 

llflflSt iMolin ari et al. 
Il996: Palla et al. 1991) . Furthermore, the detection of SiO 
maser emission confirms that the source is a late-type star, 
since SiO maser emission is extremely rare in star-forming 
regions, with only t hree (extremely luminous) sources de- 
tected to date fe.g.lE^gels fc Hesk3ll989HSnvder fc Buhll 



ll974HTJkita et alJl987|^ . Though both the 12/im and 25^m 
IRAS flux densities have good quality, their association 



Fig. 16. MSX colour-c olour plot for our SiO maser tar- 
gets, similar to Fig. 5 of iLumsden et all l|2002|) . Stars rep- 
resent AGB stars, open circles early post-AGB stars, fol- 
lowing the classification of Sevenster (2002.1 . Black body 
spectra follow the continuous line. A reddeni n g vec tor 
for — 40 and extinction law from i MathisI 11199(1 ) is 
given by the solid arrow, while another one obtained us- 
ing the extinction law of iLut j l)l999() is represented by 
the dashed arrow. Note that the two vectors point in the 
same direction. Our SiO targets have colou rs similar to 
the ob jects with 9.7 um si l icate emission in iKwok et alJ 
l|l997|) and lLumsden et all l)2002() . Object #99 is outside 
of the plotted region {A - E = 3.01 mag, Fe/Fa = 2.4; 
C -2.99 mag, Fd/Fc = 0.04) (see also text). 

is unreliable. In fact, the corresponding MSX source has 
measurements in the A, C and D bands consistent with 
the 12 /ini IRAS detection, but is not detected at 21 /im. 

7. MSX colour-colour diagrams 

ISevensteil ()2002j) analysed the mid-infrared properties of 
her OH/IR sample using IRAS and MSX data. Studying a 
possible correspondence between regions in the two-colour 
IRAS diagram and the MSX A — C versus D — E plane 
she suggested that the MSX diagram can distinguish be- 
tween the AGB and the post-AGB phases. The transition 
from a blue (< 1.8) to red (> 1.8) A — C colour may cor- 
respond to a transition off the AGB to proto-planetary 
nebulae: the star had its last thermal pulse, and ceased to 
be variable. The transition from a blue (< 1.5) to a red 
(> 1.5) D — E colour indicates a later evolutionary transi- 
tion, when mass-loss starts to drop down of several order 
of magnitudes and there is the onset of the fast wind. 
Most of our SiO targets with a clear MSX counterpart 
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Fig. 17. MSX col our-colour plot as in Fig. 5 of 
iLumsden et alJ l)2002(l for the OH/IR sample of lSevensteJ 
l)200 2*l. Symbols are as in Fig. plus crosses which rep- 
resent late post- AGE stars, following the classification of 
ISevensterl tm^ 

show A-C <1.8 and D - E < 1.5, as expected for AGB 
stars, and comparable to the bulk of Sevenster's OH/IR 
sample. 

According to Sevenstcr' criteria, only two objects, #76 
and #99, which are both SiO maser emitters, are likely to 
be post-AGB stars. The odd colour of #99 is due to an 
extremely high flux density of 39 Jy measured in the C 
band, while the flux density in both the A and D bands 
is only ~ 1.7 Jy. However, the inspection of the MSX C 
band image does not confirm the presence of a such bright 
object and we conclude that the C photometry (despite 
its g ood fiag) of #99 is unreHable. Source #76, LPV 12- 
352 l)Glass et al.ll200l|) . is located on a region of extended 
emission which is associated with a star forming region 
this could have affected the mid-infrared 
colours of #76. In conclusion: the anomalous position of 
the sources #99 and #76 is probably due to the assign- 
ment of incorrect magnitudes. 

The A — C versus D — E colour diagram is useful to 
locate post-AGB stars, which have redder colours due to 
their colder envelopes. However, this diagram can hardly 
distinguish between different thickness of the envelopes of 
AGB stars. Using a s ample of IRAS so u rces with IRAS 
low re solution spectra l|Kwok et al.lll997t) , iLumsden et al.l 
l)2002|) showed that a different combination of MSX fil- 
ters can distinguish between circumstellar envelopes with 
silicate feature at 9.7 /im in emission and in absorption. 
Miras and OH/IR stars with the silicate feature in emis- 
sion are located below the black body line in the C — D 



versus A — E diagram, while OH/IR stars with silicate 
feature in absorption lie above this line. In optically thick 
envelopes, self-absorption causes a decrease of the flux in 
the C and A bands, leading to an increase in C — D. 
Figure El shows that our SiO targets are distributed like 
objects with the silicate feature in emission. We note that 
the MSX two colour plot is not corrected for reddening. 
However, the reddening correction makes the sources bluer 
in the C — D colour independently of the adopted extinc- 
tion law. Furthermore, in the C — D versus A — E diagram 
the distribution of the ISOGAL and MSX samples are 
similar. For comparison, in Fig.[T7|the MSX colour-colour 
plots of the Sevenster's OH/IR stars is shown, which are 
distributed over a wider and redder range of colours. 

8. Conclusion 

To increase the number of line-of-sight velocities measured 
toward the inner Galaxy we searched for 86 GHz SiO 
maser emission in a sample of 441 late- ty pe stars. While 
the r adio survey was described elsewhere ()Messineo et alJ 
I2OO2I) . this paper presents the infrared photometry of all 
our SiO targets as derived from the various accessible in- 
frared catalogues: DENIS, 2MASS, MSX, ISOGAL and 
IRAS catalogues. 

As described in IPaperT . we initially selected the SiO 
targets from the ISOGAL and MSX catalogues on the 
basis of their near- and mid-infrared colours, and their 15 
/^m magnitudes. We tried to select objects with colours 
typical of pulsating AGB stars with thin envelopes (Mira- 
like stars), while avoiding OH/IR stars and other sources 
with thick circumstellar envelopes. 

Our analysis of the targeted stars' multi-band pho- 
tometry showed that these selection criteria were quite 
reliable. A comparison between the DENIS and 2MASS 
data shows that most of them are variable stars, and more- 
over the correlation between the J and band brightness 
variations is similar to that found in loca l dust-enshrouded 
Mira variable stars l|01ivier et alJl200l|) . 

The IRAS [12] - [25] colours of the SiO targets con- 
firms that they populate mostly the region Ilia of the van 
der Veen and Habing classical two-colour IRAS diagram, 
which is a region of stars with moderate mass-loss rates 
and with silicate feature at 9.7 /im in emission. The distri- 
bution of the [12] — [25] colours of the SiO targets overlaps 
with the distribution of the colours of OH/IR stars, which 
however are distributed over a larger and redder range 
of colours (mostl y in re gion Ilia and Illb). F ollowing the 
work of S cvensten l)2002(l and lLumsden et all 42002), those 
properties can be translated and seen in the MSX C — D 
vs. A — E diagram. The SiO targets have a narrower C — D 
colour range and are located below the black-body line, 
differently from the thick-envelope OH/IR stars. 

The two subsamples, the MSX-selected objects and the 
ISOGAL selected objects have very similar infrared prop- 
erties but differ slightly in the average apparent magni- 
tude, the MSX sample being on average a little brighter. 
This difference is, however, smaller than the spread in 
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magnitudes of each subsample. In forthcoming papers we 
usually will combine the two subsamples into one. 
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